#include<iostream>
using namespace std;
#include<cmath>

#include"library.h"

int main(){
	int n=3;
	double *a = new double [n];
	double *b = new double [n];
	double *c = new double [n];
	double *u = new double [n];
	double *r = new double [n];
	a[0]=0; a[1]=1; a[2]=1;
	b[0]=2; b[1]=2; b[2]=2;
	c[0]=1; c[1]=1; c[2]=0;
	r[0]=3; r[1]=4; r[2]=3;
	
	TridiagLinear(n, a, b, c, u, r);
	cout<<"u: "; for(int i=0;i<n;i++)cout<<u[i]<<","; cout<<endl;
	delete [] a,b,c,u,r;
	return 0;
}
